Microscopic versus mesoscopic local density of states in one-dimensional localization 
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We calculate the probability distribution of the local density of states :/ in a disordered one- 
dimensional conductor or single-mode waveguide, attached at one end to an electron or photon 
reservoir. We show that this distribution does not display a log-normal tail for small i/, but diverges 
instead oc u~^'^. The log-normal tail appears if v is averaged over rapid oscillations on the scale 
of the wavelength. There is no such qualitative distinction between microscopic and mesoscopic 
densities of states if the levels are broadened by inelastic scattering or absorption, rather than by 
coupling to a reservoir. 



PACS numbers; 72.15.Rn, 42.25.Dd, 73.63. Nm 

Localization of wave functions by disorder can be seen 
in the fluctuations of the density of states, provided the 
system is probed on a sufficiently short length scale [|l],|| . 
The local density of states (LDOS) of electrons can be 
probed using the tunnel resistance of a point contact 0] 
or the Knight shift in nuclear magnetic resonance Wf, 
while the LDOS of photons determines the rate of spon- 
taneous emission from an atomic transition B. (In the 
photonic case one can study the effects of localization 
independently from those of interactions.) 

For each length scale 5 characteristic for the resolution 
of the probe, one can introduce a corresponding LDOS 
Vs. It is necessary that S is less than the localization 
length, in order to be able to see the effects of localiza- 
tion — the hallmark M being the appearance of loga- 
rithmically normal tails oc exp(— constant x In 1^5) in the 
probability distribution P{vs)- 

Much of our present understanding |7| of this problem 
in a wire geometry builds on the one-dimensional (ID) 
solution of Altshuler and Prigodin |0|. In the simplest 
case one has a single-mode wire which is closed at one end 
and attached at the other end to an electron reservoir. 
(The optical analogue would be a single-mode waveguide 
that can radiate into free space from one end.) In ID 
the localization length equals twice the mean free path /, 
which is assumed to be large compared to the wavelength 
A. One can then distinguish the microscopic LDOS v = 
vs for 5 <C A, and the mesoscopic LDOS v — vg for 
A <C (5 ^ /. While v oscillates rapidly on the scale of the 
wavelength, v only contains the slowly varying envelope 
of these oscillations. Altshuler and Prigodin calculated 
the distribution P{v) and surmised that P{v) would have 
the same log-normal tails. We will demonstrate that this 
is not the case for the small-j^ asymptotics. 

The calculation of Ref. [gj was based on the Berezinskii 
diagram technique B , which reconstructs the probability 
distribution from its moments. (An alternative approach 
pO| , using the method of supersymmetry, also proceeds 
via the moments.) An altogether different scattering ap- 
proach has been proposed by Gasparian, Christen, and 



Biittiker |T^, and more recently by Pustilnik [Q. We 
have pursued this approach and arrive at a relation be- 
tween V, V and reflection coefhcients. This allows a di- 
rect calculation of the distributions. We find that P(y) 
and P(y^ have the same log-normal tail for large densi- 
ties, but the asymptotics for small v and v is completely 
different. The strong fluctuations of v on the scale of 
the wavelength lead to a divergence P{v) oc v^^l"^ for 
j^ ~> 0, while the distribution of the envelope vanishes, 
P(y) -^ for z^ ^ 0. This qualitative difference between 
microscopic and mesoscopic LDOS is a feature of an open 
system. Both P(y) and P(y) vanish for small densities if 
the wire is closed at both ends and the levels are broad- 
ened by inelastic scatterers (for electrons) or absorption 
(for photons). 

We consider a ID wire and relate the microscopic and 
mesoscopic LDOS at energy E and at a point x = to 
the reflection amplitudes rpj, tl from parts of the wire 
to the right and to the left of this point. The Hamilto- 
nian is H = —{'h?/2m)d^/dx'^ + V{x) for non-interacting 
electrons. (For photons of a single polarization we would 
consider the differential operator of the scalar wave equa- 
tion.) We will put h= \ for convenience of notation. We 
start from the relation between the LDOS and the re- 
tarded Green function. 
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with r] a positive infinitesimal. We assume weak disorder 
(kl ^ 1, with k = 27r/A the wavenumber), so that we 
can expand the Green function in scattering states in a 
small interval around x = 0, 



G(a;) = CL(e-^'^" + rLe*^'")0(- 
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(The function d{x) ^ 1 ior x > and for a; < 0.) The 
coefficients cl and cr are related by the requirement that 
the Green function be continuous at x = 0, 



cl(1 + r-L) = cr(1 + rR). 



(4) 



Substitution of Eq. (||) into Eq. (||) gives a second relation 
between cl and cr, from which we deduce that 



G(0) 



(l + rL)(l + rR) 
iv{l - TRrL) 



(5) 



with V the velocity. Using Eq. (nh we arrive at the key 
relation between the microscopic LDOS and the reflection 
coefficients, 

i^ = (ttv)-^ Re (f + rL)(l - rRr^y^l + tr). (6) 

In order to perform the local spatial average that gives 
the mesoscopic LDOS D, we use that the reflection co- 
efficients oscillate on the scale of the wavelength. If we 
shift xq slightly away from the origin to a point x' , one 



has tl -^ g^«a; ^^ g^^^ j,^ 



'2ih. 



^ tr. The product 



rRTh, however, does not display these oscillations — only 
this combination should be retained. Hence 



j> = (ttw) ^Re(I + rRrL) (f - TRrL) ^ 



(7) 



In what follows we will measure v and i> in units of vq — 
(7rw)~^, which is the macroscopic density of states and 
the ensemble average of i^, />. 

Let us now demonstrate the power of the two simple 
relations (||) and (0). We take the wire open at the left 
end and study the density at a distance L from this open- 
ing. At the right end the wire is assumed to be closed, 
giving rise to a reflection coefficient tr = exp(i0R) with 
uniformly distributed phase 0r in the interval (0, 27r). 
The reflection coefficient tl — vi?exp(i0L) is parame- 
terized through the uniformly distributed phase (f>h and 
the reflection probability R in the interval (0,1). The 
ratio 



u= (l + i?)(l-i?)-^ 
has the probability distribution ||l3| 
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with s = L/l and I the mean free path for backscattering. 
The mesoscopic LDOS (0) can be written in terms of the 
variables u and (p — (pi, + (fiji, 



D ^ iu — y V? — 1 cos ) 



(10) 



Averaging flrst over </> wc find 



Popenii>) ^ '^^^ f du -^^, a =i(l> + !>-!). (11) 
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The subsequent integration with Eq. (m yields 
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The distribution function (n3) is the celebrated result 
of Altshulcr and Prigodin B . It displays log-normal tails 
for both large and small values of D. Indeed, the two tails 
are linked by the functional relation M 



P{l/v) = vT{v) 



(13) 



This relation follows directly from Eq. (|11|) and hence re- 
quires only a uniformly distributed phase </>, regardless of 
the distribution function p{u) of the reflection probabil- 
ity. As we will now show, such a relation does not hold in 
general for the microscopic LDOS v, and the asymptotics 
of its distribution function for small and large values of 
V can be entirely different. 

The calculation is facilitated by the fact that v is re- 
lated to V by 



i/ = 2z>cos^(</«r/2) if|rR| = l. 



(14) 



Moreover, v is statistically independent of 0r because 
the latter enters v only in combination with (/)l, which 
itself is uniformly distributed. The distribution of the 
microscopic LDOS hence follows directly from Eq. (llT 
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where we substituted t — cos^(0r/2). The asymptotic 
behavior is 



Popcn(^) - 21/2^ "" 



exp(3s/4) _ _i/2 



J^< e" 



(16a) 



-* open v^j 



2V2exp(~s/4) 3/2 
si/27r3/2 



e-'<^v<^e', (16b) 



_ exp[-./4-ln^(./2)/4s] 
°^°"^^" 7r.3/2ini/2(,/2) ' "»^- ^'^'^ 

In the second and third region this is similar to the be- 
havior of -Popcn(^') in Eq. (12). In the region of the 



smallest densities, however, Popen('^) is not log- normal 
like Popon(i') but diverges oc z/~^/^. 

The different tails arise from two qualitatively differ- 
ent mechanisms that produce small values of v and v. 
For the mesoscopic LDOS this requires remoteness of E 
from the eigenvalues of wave functions localized within a 
localization length around xq. As a consequence, P{v) 
is intimately linked to the distribution function of reso- 
nance widths IQ]. Small values of the microscopic LDOS 
V are attained at nodes of the wave function which solves 
the wave equation with open boundary conditions, inde- 
pendent of the energy. The nodes are completely deter- 
mined by the small-scale structure of the wave function. 
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FIG. 1. Distributions of the microscopic local density of 
states (LDOS) v and the mesoscopic LDOS v for the open 
wire at a distance L = 21 from the opening. (Both are mea- 
sured in units of their mean uo = {ttv)~^.) Sohd curves are 
given by Eqs. (hd) and ([L5|). The data points result from a nu- 
merical simulation for a wire of length 10 I with no adjustable 
parameter. The inset shows the geometry of the open wire 
(not to scale). 



which is a real standing wave ex cos {kx + a) with ran- 
dom phase a (t) . [Wc recognize the square of this wave 
amplitude in Eq. ([14]).] The resulting v~^^'^ divergence 
of the probability distribution has the same origin as in 
the Porter-Thomas distribution for chaotic wave func- 
tions ll]. 

The two distributions for the open wire are plotted in 
Fig. n^, together with the result of a numerical simulation 
in which the Green function inside the wire is calculated 
recursively p5| | . The comparison of theory and numer- 
ics is free of any adjustable parameter — the velocity was 
taken from the dispersion relation, and the mean free 
path was obtained from the disorder strength within the 
Born approximation. 

We now show that this qualitative difference between 
the microscopic and mesoscopic LDOS is absent in a 
closed wire. If the wire is decoupled from the reservoir 
we need another source of level broadening to regularize 
the delta functions in the LDOS. Following Ref. S], we 
will retain a finite imaginary part r\ of the energy, cor- 
responding to spatially uniform absorption (for photons) 
or inelastic scattering (for electrons), with rate Ij]. Eqs. 
(^ and (^ still hold provided i] <$i E. The reflection 
coefficients can be written as 



''R.L = \/RR,he''' 



(17) 



where (/)r and (/)l are uniformly distributed phases and 
-Rr, -Rl are independent reflection probabilities. In an 
infinitely long wire they have the same distribution p6[ 
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■exp[-Lj{l- R)-^], u; = 'iT]l/v. (18) 
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FIG. 2. Same as in Fig. hi but for the closed wire with 
dimensionless absorption rate oj = 1/6. Solid curves are given 
by Eqs. (Elf) and (B3). The data points result from a numerical 
simulation for a wire of length 55 1, with the LDOS computed 
halfway in the wire. 



After elimination of the phases the distribution of the 
mesoscopic LDOS takes again the form (11), where u now 
stands for the combination 



M=(l + i?Ri?L)(l-i?Ri?L)"'. 

Eq. ( |l8| ) implies for u the distribution 



(19) 



(20) 



The resulting distribution function of the mesoscopic 
LDOS is 



Pc 



closed 



(^) 



i>-3/2 r du 



j{u-l} 



7r\/2 Ja \/u - a, 
uKo (lu^/u^ - l\+^/u^ - IKi (Loy/u"^ - l\ , (21) 

with a defined in Eq. (O) . It vanishes for small densities 



^closed (i^) = 2 ^'^ujv 2exp(-w/£>), v<^u. (22) 

This should be compared with the known distribution |^ 

/'cioscd(i^) = Y ^^"'^' exp [u - \u:(v + v-^)\ (23) 

of the microscopic LDOS. In contrast to the open wire, 
both distributions vanish for i/, i/ — > 0. This is illustrated 
in Fig. 0, which compares the analytical predictions to 
numerical data obtained by diagonalization of a Hamil- 
tonian. The comparison is again free of any adjustable 
parameter. 

We note in passing that the asymptotic behavior (p2 
differs from the asymptotic behavior 



Pcioscd(j>) + liTTioY^^D-^/^ exp(-^2^/16i>), (24) 

given in Ref. [g| for oj <C 1. There the distribution func- 
tion was reconstructed from the leading asymptotics of 
the moments Uni;j_+o(t'") = w^~"n!/(2n— 1). This would 
be a valid procedure if the distribution would depend only 
on the product Lui> in the limit uj —^ 0, which it does not. 
The subleading terms of the moments have to be included 
for i> < Lu. Indeed, our distribution function has the same 
leading asymptotics of the moments, but has a different 
functional form. This illustrates the potential pitfalls of 
the restoration procedure which are circumvented by our 
direct method. 

In conclusion, we have given exact results for the distri- 
butions of the local densities of states in one-dimensional 
localization, contrasting the microscopic length scale (be- 
low the wavelength) and mesoscopic length scale (be- 
tween the wavelength and the mean free path) . Contrary 
to expectations in the literature, the log-normal asymp- 
totics at small densities applies only to the mesoscopic 
LDOS !/, while the distribution of the microscopic LDOS 
v diverges ex v~^l'^ for i/ ^- 0. This is of physical signif- 
icance because many of the local probes act on atomic 
degrees of freedom and hence measure v rather than v. 
The strong length scale dependence of the LDOS disap- 
pears if the electrons (or photons) are scattered inelasti- 
cally (or absorbed) before reaching the reservoir. Both 
Piv) and P(y) then have an exponential cutoff at small 
densities. 

It is an interesting open problem whether the qualita- 
tive distinction between v and v in an open wire carries 
over to the quasi-one dimensional geometry with iV > 1 
modes. An analytic theory could build on the multi- 
channel generalization of Eq. (n) , 
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ReTrM(l + fi^){l - fnfLy\l + tr). 



(25) 



Now tl and f^ are N x N reflection matrices and the 
matrix M„m = 2(7rA)"i(w„i;m)~i/2 sin(g„ •ro)sin(g„ -ro) 
contains the weights of the TV scattering states with 
transversal momentum (f„ and longitudinal velocity u„ 
at the transversal position tq on the cross section of the 
wire (area A). 

Another promising direction for future research is to 
study what happens to the LDOS if the wire is coupled 
to a superconductor rather than to a normal electron 
reservoir [l7| ]. The convenient expressions for v and D 
in terms of the reflection matrices from two independent 
parts of the wire, derived in this paper, can be directly 
generalized to include Andreev reflection at the interface. 
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